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Abstract 

We present the first numerical QCD bound state calculation based on a renormalization 
group-improved light-front Hamiltonian formalism. The QCD Hamiltonian is determined 
to second order in the coupling, and it includes two-body confining interactions. We make 
a momentum expansion, obtaining an equal-time- like Schrodinger equation. This is solved 
for quark-antiquark constituent states, and we obtain a set of self-consistent parameters 
by fitting B meson spectra. 
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1 Introduction. 

Recently, a new approach to renormalization and solving for QCD bound states has been 
advocated [||, |2|, ^, Q. The goal is to build a bridge between QCD and a constituent quark 
model (CQM). It has been argued that it is convenient to use a light-front formulation 
of the theory, because on the light-front it is possible to make the vacuum trivial simply 
by implementing a small longitudinal cutoff. As a result all partons in a hadronic state 
are connected to the hadron, instead of being disconnected excitations in a complicated 
medium. The price to pay is a considerably more complicated renormalization problem. 

In this paper we briefly describe this new approach, concentrating on aspects, necessary 
to appreciate the simple calculation we present. The calculations of the heavy meson 
spectra, for systems containing one heavy and one light quark, is intended to elucidate 
the approach and provide qualitative tests of the leading terms revealed in the effective 
Hamiltonian by the renormalization group. Only the first step in this calculation is taken 
here, and later work will focus on the spin-dependent structure. We refer the reader to 
literature for many details |jl], |2|, ^, ^ . 

The new approach consists of two steps. The first step is renormalization. The second 
step is a bound state calculation. 

The aim of the first step is to find an effective renormalized Hamiltonian at hadronic 
energy scales starting with a Hamiltonian which is consistent with perturbative QCD at 
high energy scales. A natural starting point is the canonical light-front QCD Hamilto- 
nian, although it cannot be complete — there must be other operators, which cannot be 
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determined from the perturbative behaviour of the theory. Glazek and Wilson de- 
signed a similarity transformation to lower the cutoff scale, which is tailored to make the 
Hamiltonian look more like a CQM Hamiltonian (see, for example, and references in 
Q). The cutoff violates manifest gauge invariance and Lorentz covariance, and thus these 
symmetries are no longer a guide to what operators are allowed in the Hamiltonian. 

Similarity transformations can be designed to bring the Hamiltonian toward a band- 
diagonal form by eliminating matrix elements between states which differ drastically in 
light-front energy. Effects of couplings that are removed have to be put directly into the 
Hamiltonian as new effective interactions. One important consequence is that two-body 
potentials are generated. In fact, the similarity transformation generates a logarithmic 
confining potential already at order ||2[0. This part of the calculation is done perturba- 
tively. If the similarity transformation can be done analytically, it is easy to use a powerful 
method called coupling coherence ^ to determine all new terms. We will illustrate this 
procedure, which is straightforward to second order, in the following sections. 

The second step is the bound state calculation itself. The effective Hamiltonian is 
divided into Hq, a part which is solved nonperturbatively, and V, the difference between 
the original Hamiltonian and Hq. The effects of V are to be computed using bound 
state perturbation theory. The criteria for choosing Hq are that it approximates the 
physics relevant for hadronic bound states as closely as possible (we take a hint from 
the CQM and include constituent masses and two-body potentials) and yet it must be 
manageable. We emphasize that the approach is tailored within limits to take advantage of 
the successful phenomenology, but it does not stop there. We can systematically improve 
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the calculations, both by computing corrections to the Hamiltonian and higher higher 
terms in bound state perturbation theory. For example, any terms added by hand (e.g. 
constituent masses) can be added in such a way that at the physical value of the coupling 
they reduce to terms in the effective QCD Hamiltonian (for details see 0, Q). 

In this paper we present one of the simplest possible calculation of QCD bound states 
based on the new approach. In the first step we find the effective Hamiltonian to order 

using a similarity transformation and coupling coherence, and in the second step we 
solve for the lowest lying qq color singlet states with arbitrary but nonzero masses in the 
nonrelativistic limit. These approximations are severe, but we will see that the qualitative 
results are good. We will not explicitly show operators that have zero expectation value 
in the qq color singlet. 

We wish to derive an effective Hamiltonian that acts at the hadronic scale by lowering 
the similarity cutoff perturbatively as low as possible, hopefully down to the scale at which 
the bound state is well approximated by its two particle component, qq. This may not be 
possible for all systems. The coupling may become too large for perturbation theory to 
be reliable before the higher Fock components are eliminated. We know, however, that 
it is possible in QED, and we believe that the success of the CQM suggests the same for 
QCD. The most favorable systems as far as the qq and nonrelativistic approximations are 
concerned are heavy quarkonia. However heavy mesons containing a light quark provide 
a better tool for testing our approach qualitatively, because they have fine structure that 
allows us to separately test spin and orbital angular momentum-dependent operators 
in the Hamiltonian, in the regime which is dominated by the confining interaction. This 
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means that the scahng of momenta and energies is determined by the confining interaction. 
As the Coulomb part of the potential becomes more important, the simple scaling analysis 
breaks down, until the Coulomb potential dominates over the confining potential. Then 
the momenta and energies scale as in QED. We expect heavy quarkonia to be in the 
mixed regime, so the simple scaling analysis is probably not going to be reliable. Study 
of the spin and orbital angular momentum-dependent operators, which are generated by 
the similarity transformation to second order in the coupling, will be done in later work. 

In heavy-light systems there is only one heavy quark, but we do not yet know whether 
the light degrees of freedom can be approximated by just one constituent in our approach. 
Further, as we show below, heavy-light mesons are qualitatively different from heavy 
quarkonia, but there are similarities with other mesonic systems: strange mesons and 
isospin 1 light mesons. So some of what we learn from heavy mesons may be generalized 
to light mesons. We will choose B mesons to check whether we can fit spectra with 
reasonable parameters. 

After bringing the cutoff down to the hadronic scale, we do bound state calculations. 
We have to decide which terms will be treated as the dominant interactions and put in 
the unperturbed Hamiltonian, and which ones will be treated as perturbations. 

In order to gain qualitative insight it is useful to study the nonrelativistic limit and to 
rewrite the bound state equation in position space. It turns out that in the nonrelativistic 
limit light-front dynamics naturally reduce to equal-time dynamics, which implies that 
angular momenta become kinematically defined. It is possible to transform light-front 
coordinates to equal-time coordinates by a specific change of variables |0] without tak- 
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ing the nonrelativistic limit. However, the point we want to make here is that in the 
nonrelativistic hmit equal-time dynamics arise naturally. This will become clearer below. 

Nonrelativistic reduction can be justified at best only for the lowest lying states. We 
may need to do a series of bound state calculations with small coupling, and then extrapo- 
late to the physical value of the strong coupling if this is large However, we are not yet 
at the stage where we can carry out the strategy with confidence that we have complete 
control over the bound state perturbation theory. The similarity transformation generates 
effective operators. The primary motivation of this work is to initiate the study of these 
operators. Further, the issue of chiral symmetry breaking terms, that is, what operators 
have to be added to the Hamiltonian by hand to restore the effects of zero modes removed 
by a small longitudinal cutoff, is not resolved yet. Nonrelativistic reduction provides a 
framework in which these and other questions can be studied. 

The main questions answered in this calculation are: 

1. Two-body potentials generated by the similarity transformation include the Coulomb 
potential and a logarithmic confining potential already at order g^. Does this effec- 
tive Hamiltonian contain enough structure to provide a starting point for studying 
bound states in QCD? For the spin-independent part of the effective Hamiltonian, 
which we study here, the answer is: Yes. 

2. The similarity transformation generates new effective interactions. Is there a simple 
way to see how these new operators affect spectra? Yes. 

3. In order to answer the previous questions, one has to make some simplifications of 
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the original Hamiltonian, including but not restricted to nonrelativistic reduction 
(for details see section 3). We do not focus on quantitative analysis, but we want to 
ask whether the calculation is at all reasonable; that is, can we fit any data with a 
set of reasonable parameters? Yes. 

The paper is organized as follows. In the second section we briefly review the general 
strategy, sketching the similarity transformation and coupling coherence. The third sec- 
tion gives our calculation - qq to order . First, we find the effective Hamiltonian using 
the similarity transformation and coupling coherence. Then we split the Hamiltonian into 
a part which is treated nonperturbatively in the bound state calculation and a part which 
is treated perturbatively. In this paper, we solve the leading order problem, and show 
that it reduces to a simple Schrodinger equation in the nonrelativistic limit. We show a 
simple way to qualitatively analyze the physics behind this Schrodinger equation. Then 
we present the numerical results. The last section contains our summary and conclusions. 

2 Two steps to solving QCD - general strategy 

In this section we briefly review the general strategy, first outlined in ref. Q. We start 
with the canonical light-front QCD Hamiltonian in light-cone gauge, = 0. We will 
not explicitly show terms that are not important for the specific calculations presented 
in the next section. For a detailed discussion of the light-front Hamiltonian see |ll| 
and references therein. Ignoring purely gluonic terms that do not affect the effective qq 
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Hamiltonian until fourth order in g, 



H = Hiree + Vi+V2 ; 



(1) 



where .fffree is the free hght-front Hamiltonian: 
H 



free 2^ / (27r)32p+ p+ 



= E 



{b\bf+d\df) . 



(2) 



There is a sum over flavors, is a quark mass, bf, df are quark and antiquark annihilation 

operators; 

Vi= g j dx'd'^x±i>4.tp (3) 

contains the standard order g quark-gluon coupling. Here ip and A^^ = J2a ^a'^"' are free 
light-front fields: 



\ 



and 



= = 0, A- = ■ A^, A^ 



The constrained fields, V- and A , are replaced by functions of the physical degrees of 
freedom resulting in new terms in the canonical Hamiltonian, among which 



V2 = -2^1 dx-d^x±{^iT'^^+) 



(4) 



is the instantaneous gluon exchange between two fermions. 



We regulate the Hamiltonian with cutoffs on the change in free energy at each inter- 
action vertex and a cutoff on longitudinal momentum fractions that is taken to zero at 



the end of the calculation. Then we use a similarity transformation The similarity trans- 
formations form a renormalization group. Repeated transformations generate a sequence 
of Hamiltonians with decreasing cutoff. A Hamiltonian in the sequence is related to the 
previous one by: 



U is a unitary matrix and can be written as U = e*^, where R is hermitian and has an 
expansion in powers of the nondiagonal part of the Hamiltonian. We have chosen n so 
that it decreases as the cutoff decreases. 

The similarity transformation is designed to bring the cutoff Hamiltonian to a band 
diagonal form, while avoiding problems of small energy denominators in perturbative 
expansions |5|. In particular, the transformed Hamiltonian is required to be band diagonal 
relative to the new scale. This means that the matrix elements between states that differ 
in energy by more than the new cutoff must be zero for the simple step function cutoffs 
we employ here. This requirement has implications for the matrix elements of R. Given 
R, one can find the transformed Hamiltonian. 

The initial cutoff destroys manifest symmetries and one of the criteria for the renor- 
malized Hamiltonian is that it restores the symmetries, albeit not necessarily in manifest 
form. If the similarity transformation can be done analytically, it is possible to use cou- 
pling coherence ^ to completely fix the renormalized Hamiltonian without explicit 
reference to symmetries. The basic idea of coupling coherence is that in the Hamiltonian 
restricted by symmetries, the strengths of all operators are not independent but depend 




(5) 



9 



only on a finite number of independent canonical parameters; so that under a transfor- 
mation, the Hamiltonian reproduces itself in form exactly, apart from the change of the 
explicit cutoff and the running of those few independent couplings. All dependence on 
the cutoff is absorbed into the independent running couplings. Once one obtains a Hamil- 
tonian that reproduces itself as the cutoff is lowered, any initial cutoff can be sent to 
infinity. 

In order to use coupling coherence we need to study how the Hamiltonian changes 
when the cutoff changes. Let = //free + v, where i/free is a free Hamiltonian and 
V = Vi + V2 'is cut off so that 

{<l)M(t>j) = , (6) 

if \EQi — E'ojl > p^; where i/free I = ^Oil^i)- The similarity cutoff, with the 
dimension of light-front energy, consists of A^, which carries a dimension of transverse 
mass squared, and an arbitrary longitudinal momentum reference scale . If this cutoff 
is lowered to by the similarity transformation, the new Hamiltonian matrix elements 
to are 0: 



HA„-iab = {(t>a\ Hfj-ce + V 



T,k VakVkb 



^2 .2 



e{\A^k\~^)e{\Aat:\-\A,,\) e(|Ai,fc|-^)e(|A,fc|-|A,fc|) 

Eok — Eoa Em^—Eojj 



(7) 



A2_ 

where Aij = E^i — Eqj and l-Eoa — -E'obI < ; and there are implicit cutoffs in this 
expression because the matrix elements of v have already been cut off so that Vij = if 
\Eoi — Eoj\ > 

To this order, a coupling coherent Hamiltonian reproduces itself, with the only change 
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being A„ A„_i. At the third order one begins to see the quark-gluon couphng run. If 
the Hamiltonian reproduces itself, the index n becomes irrelevant. The solution is found 
by noting that we need the partial sum above to be added to an interaction in v that is 
expressed as a sum, so that the transformation merely changes the limits on the sum in 
a simple fashion (for details see P]). There are two possibilities. The first is: 



H, 



ab 



h\ho + v\ 



- Efc VakVkb 



and the second is: 



_E± 

Eok-Eoa 



+ 



(8) 



Hab = {(f>a\ho + v\(f)b) 
+ T,k VakVkb 



g( ^-|A_,fc|)e(|A„,|-|Aft,|) ^ e(^-|A,fc|)e(|Aftfe|-|A„fc|) 



Eok — Eoa 



Eok.—En 



(9) 



Note that v in these expressions is the same as that above only to first order. The 
coupling coherent interaction in H is written as a power series in v which reproduces 
itself under the transformation, except the cutoff changes. In higher orders the canonical 
couplings also run. To decide which of these equations to use one must in principle go to 
higher orders, but in practice it is usually obvious which choice is correct. For example, 
the first solution provides effective two-body interactions from one-gluon exchange; while 
the second provides the relevant part of the quark self- energy. If chosen otherwise, the 
new terms would make the effective Hamiltonian divergent. 



At the end of this first step, the Hamiltonian is renormalized and the scale relative 
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to which it is band diagonal is a hadronic scale. [| The effective Hamiltonian contains 
complicated potentials, which result from eliminating the coupling between high and low 
energy states. It still contains emission and absorption interactions, but these no longer 
mix states of high and low energies. 

In the second step, the Hamiltonian which is now given as an expansions in g is 
regrouped from the point of view of what is important in the bound state. Based on 
the success of the CQM, we believe that emission and absorption processes, which would 
mix different Fock states, are suppressed with respect to interactions that do not change 
particle number for low-lying states and sufficiently low cutoff. Thus the particle number 
changing interactions will be treated perturbatively using bound state perturbation theory. 
This is clearly justified if the gluons are massive but even for massless gluons one can 
argue that the interacting gluon effectively acquires mass related to the confining scale 
Q . It is therefore plausible to assume that adding an extra constituent is suppressed even 
in this case. 

Next, we use constituent masses. It is possible to add the constituent mass at zeroth 
order and subtract it at a higher order in bound state perturbation theory, as outlined 
in Q. This issue starts to be important when the effective Hamiltonian is calculated to 
higher orders and one tries to see that the approximations in the leading order bound 
state calculation lead to convergent results. 

In the calculation presented here we also use a nonrelativistic reduction and we choose 
A specific value of the "hadronic scale" is to be specified by fitting spectra. 
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a rotationally symmetric Hq. We want to minimize higher order corrections, but we do 
not yet know what restores the rotational symmetry (e.g., adding extra gluons or terms 
of higher order in g). The answer to this question may change our choice of Hq in future 
calculations. At least for now, we choose the confining potential such that it does not yield 
any corrections in first order bound state perturbation theory. The corrections which come 
from the rotationally noninvariant terms are not suppressed by powers of the coupling, 
despite the fact that they enter at the second order of bound state perturbation theory. 
This is the same order at which one has to include emission and absorption processes, 
unless those are suppressed nonperturbatively (e.g. if the gluons are massive). Therefore, 
it is not clear whether it is meaningful to consider the corrections due to rotationally 
noninvariant terms without including qqg. Nevertheless, we evaluated these corrections 
for the ground state and they are of order a few percent in a region where the confining 
potential dominates over the Coulomb potential. 



3 The simple QCD calculation - qq to order 

In this section we discuss mesons, i.e. color singlet QCD bound states whose valence 
constituents are a quark and an antiquark. The masses of the constituents are arbitrary 
but nonzero. We expect that the qualitative aspects of the study are relevant to all qq 
systems with a possible exception of light isospin zero mesons^. We fit B mesons. 

To order g^, the similarity transformation is represented by a few diagrams, as shown 
^For light mesons with 1=0 we expect that operators of 0{g*) will play an important role. 
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fermion in a high energy state 
gluon in a high energy state 



Figure 1: Diagrams representing similarity transformation to order for qq. The top line 
represents new effective one-body operators, the bottom line are two-body operators. We use 
untypical lines for gluons and fermions to emphasize that these represent new operators, not 
Feynman diagrams. 

in figure 1, so it is possible to find the effective Hamiltonian analytically. Let us work out 
in detail one of the operators, and then list the results for remaining ones. 

Let us consider matrix elements involving one gluon exchange between a quark and an 
antiquark (see the bottom two diagrams in figure 1). If the Hamiltonian is band-diagonal 

14 



relative to a scale then only those matrix elements in which the energy difference is 
less than this value are nonzero: 

-9A,u{P2,cr2h^u{pi,ai)v(k2, X2h''v{ki, Xi){TaTb) 



0(J|-|I?i|M|Di|-|D2|) , e{^-\D2\)e{\D2\-\D^\) 



+ 



D2 



H-g+) 



n , ^ / S(^~\-Dim\Di\-\D2\) , e{^-\-D2mD2\-\Di\) 
-^Mi'V D I -Di + -D2 



(10) 



where pi, ki are light-front three- momenta carried by the constituents; cjj, Aj are their 



light-front helicities; u{p,a), v{k,X) are their spinors |12]; index i = 1,2 refers to the 



initial and final states, respectively; D^y{q) = ^r/^r/,, + 4- [VtiQ u + Vuq\] - 5, 



IS 



the gluon propagator in light-front gauge [1^], r/^ = (0, = 1,0,0); q = pi — p2 is the 



exchanged momentum, (7^ = ^^; Di, D2 are energy denominators: Di = Pi — P2 — OT 
and D2 = k2 — k^ — q~ . 

It is convenient to use Jacobi momenta. Setting the total transverse momentum to be 
zero, the momenta of the constituents are: 



pi = XiP^ , pf- = Kf- ; 
kf = {l — Xi)P^ , k^ = —nf ■ 

Let the mass of the constituent with momentum p be rua and the mass of the other 
constituent be nih. The denominators in terms of Jacobi momenta are: 



P+ I Xl X2 Xl — X2 

1 / /t2"^ + ml Ki'^ + ml (ki - /?2")^ 



P+ \ 1 — X2 1 — Xl 



Xl - X2 



(11) 



15 



When the scale is lowered to -p^ by the similarity transformation, all matrix elements 
in which the energy jump is larger than the new cutoff are zeroed. The effects of couplings 
which are removed have to be put directly in the Hamiltonian as new effective interactions. 
In this case, the new effective interactions according to (8) are: 



e{q+) 



-gl^_^u{p2,a2)Yu{pi,ai)v{k2A2)Yv{ki,\i){Tan) 
n , ^ ( 9{\Di\-^)e{m-\D2\) , e{\D2\-^)e{\D2\-\Di\) 



+ ^zy D^,{-q) I 2±_ + E±_ 



(12) 



This will repeat as the cutoff is lowered. Once the interaction reproduces itself in form 
the initial cutoff can be sent to infinity and we can lower the cutoff to the scale of interest. 
However, as the cutoff decreases, the coupling increases and at some point it becomes 
invalid to use perturbation theory to further lower the cutoff scale. The interaction in the 
effective Hamiltonian at the hadronic scale A is thus: 



-9au{p2, (72)l^u(j)i,ai)v{k2,\2)Yv{ki,\i){TaTb) 

1 n ( J e{\D^\-^)e{\D^\-\D2\) , 0{\D2\-^)9{\D2\-\D^\) 

^^MD Di + Di 



(13) 



Here we summed the two terms corresponding to the two time-ordered diagrams in figure 1. 

Similarly, one can find effective one-body operators (self-energies) (see the top two 
diagrams in figure 1): 



2nP+ 1 r+ 



log 



eV+ 



+ 2^A2log 
3P+ 



V 



.P+ 



1 m 



V 



2V+ 2 xgiA2 + m? 
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X 



2 



(14) 



Here m and x stands for either rua and Xa, or rrib and Xh- The first term is infrared 
divergent (e is an infinitesimal cutoff on a longitudinal momentum taken to zero at the 
end of calculation) but it exactly cancels with the infrared divergence in the effective 



By finding these counterterms we have completed the first step of the calculation. Let 
us summarize the effective Hamiltonian: 



where .fffree is the kinetic energy; Vi is 0{g) emission and absorption with nonzero matrix 
elements only between states with energy difference smaller than the hadronic scale 
V2 is 0{g'^) instantaneous interaction (with no cutoff) and V2 eff includes the effective 
interactions, also 0{g^), given in previous formulae. 

This brings us to the second step: we have to regroup and approximate the Hamil- 
tonian for the purpose of bound state calculations. As mentioned above, the emission 
and absorption are not included in Hq. We include kinetic energy, instantaneous fermion 
exchange, self-energies and the most infrared divergent piece of the effective interaction 
arising from one gluon exchange. Even this is still quite complicated as a starting point 
to gain intuition, so we consider a nonrelativistic limit of this Hamiltonian. As before, we 
derive in detail results for a specific operator and list the results for the remaining ones. 

Let us consider the most infrared divergent piece of the effective interaction arising 



two-body operator (13), if the qq pair is in a color singlet 0|. 



HeS — -fffree + ^1 + ^2 + ^2 eff , 



(15) 
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from one gluon exchange in a color singlet: 



1 ( e{\D,\-^)e{\D,\-\D2\) e{\D2\-^)e{\D2\-\Di\) 



(16) 



Let us denote M^b = iria+mh, introduce a new variable for the longitudinal momentum 
fraction rj such that Xa = ^ — n, xt, = 1 — Xa and then make an expansion in powers 
of 77. 

To the lowest order in momenta, both energy denominators reduce to 

i'.^i'.^-^(Mita-^)+<|^). (17) 

Here {nj- — )^ = Q^"^- further identify = M'^^{rii — r/2)^, then the expression 

(16) reduces to: 

-.-L2 



4^1^X1(1 - XI)X2{1 - X2)Mli,\^{l - Obelow) , (18) 

where 

_ / A2 Mab<f \ 



^below — 



Prom now on, we drop the omnipresent y/xi{l — a;i)x2(l — 0:2), because it cancels exactly 
with the same factor in the definition of the wavefunction (see eqn. (30) in the next 
section) . 

Adding this to a canonical term of order which has a similar structure, namely the 
instantaneous interaction: 

1 M2 

-AglCp- =-4glCF^ , (19) 
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leads to the following interaction: 



7tl r,2pfi 



Shelow ■ (20) 



The first term is the Coulomb potential. The scale dependent part of the interaction (i.e. 
the second term) leads to a logarithmic confining potential 0] . 

Similarly, one can show that the kinetic energy in the nonrelativistic limit reduces to 

-^' + --1^ + < _ 2A4,^ , (2f) 
Xa Xb 2m 

where m is the reduced mass; and that the self-energy produces only a constant shift: 

aCpMabC 



2 log 



log ( ) + fl + ^ Vog f -^^) - - 



(22) 



where 



A2 p^ 

£ = 



r+Mab 

carries dimension of mass and in the nonrelativistic limit replaces the light-front cutoff 
Note that we did not assume anything about the relation between masses and the 
cutoff C The assumption that momenta are small in comparison to masses is equivalent 
to assuming that the wave function is peaked at small momenta, which typically requires 
that g\ is sufficiently small. 

It is more intuitive to work in position space. The Fourier transform of the potential is: 

r 

_ 2MabCFaC fi l-t_ f IjoiVt^Cr^) + J2iVt~fiCr^)] - l] 

TT Jo t ^ I- J J 

^2_AUC^ r ^eit -t^-w) {cositCrMV^Cr^) " 1} (23) 
vr Jo t"^ + w 
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Figure 2: Violation of rotational symmetry in the confining potential. At small distances, the 
violation of rotational symmetry is small. At any fixed value of 7^, the confining potential is 
maximal when the quarks are separated in purely transverse direction. It is minimal when the 
separation between the quarks is in purely longitudinal direction. We also show the strength 
of the potential averaged over the angle 9. 



where Vz is the z-component of the separation between the quarks and r± = |rj_| is the 
transverse separation between the quarks. Integration variables t and w are dimensionless. 
The first term in the expression is the Coulomb potential, the rest is the confining potential 
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normalized to zero at the origin. As mentioned above, for the color singlet states, the 
infrared divergence in this two-body operator precisely cancels with the infrared divergence 
in the self-energies, making the confining potential finite at the origin. Any finite terms 
required to make the confining potential vanish at the origin are subtracted from the 
self-energies. 

The confining potential is not rotationally invariant, because the cutoff violates ro- 
tational symmetry. Also, recall that rotations on the light-front are not kinematic, so 
as long as the gluon emission and absorption is allowed, complete rotational invariance 
requires states containing arbitrarily large numbers of gluons. The renormalized Hamilto- 
nian restores rotational symmetry, but only to 0{g^), while in the bound state calculation 
the confining interaction is treated to all orders. 

The expansion of the confining potential in Legendre polynomials has only even terms: 




k 



(24) 



with 



Vo = 21og7^-2C^(7^)+4^-2^i^^ 



+ 2 



sin 7^ 



-5 + 27 , 



(25) 



5 55^(7^) _ 10 5Si{n) BcosU 
15 bcosTZ bsinlZ 5sinTl 20 



20cos7^ 



(26) 



7^3 27^3 7^3 27^4 



and 



V2k 




TT 



V{f^P2k{x) 
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where Tl = Cr\ Ci{x), Si{x) are cosine and sine integrals, respectively; and 7 is the Euler 
constant. 

When the separation between the quarks is purely in the z-direction, i.e. r± = 0, the 
potential has its minimum value with respect to the angle 9. Integrals in (21) can be done 
analytically leading to: 

2MabCFaC 



TT 



log(7^) - C^(7^) + + -j^ - (2 - 7) 



(27) 



The potential is maximal in the purely transverse direction (i.e. = 0). For large 
separations, the difference between the potential in purely transverse and purely longitu- 
dinal directions is a factor of two. At very short distances, even the confining potential 
is rotationally invariant (see figure 2). This is of no help though, because at very short 
distances, the Coulomb part of the potential dominates. 

For the bound state calculation, we want to choose a rotationally invariant confining 
interaction for Hq. This is partly motivated by phenomenology and partly by our desire to 
find a tool to qualitatively analyse the effective operators in the Hamiltonian. Restoration 
of rotational symmetry in Hq provides such a tool in the form of a simple scaling analysis. 

It is not clear how we should choose the leading rotationally invariant interaction at 
order g^, because we do not know yet how rotational symmetry is restored. 

It is clear from (24) that it is convenient to use the first term in the expansion in 
Legendre polynomials. It is the only choice that does not lead to any corrections in the 
first order of bound state perturbation theory. The corrections which come from the 
rotationally noninvariant terms enter at second order in bound state perturbation theory; 
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that is, at the same order as the emission and absorption processes. We have to mention 
that these corrections are not suppressed by powers of the coupling, and that the problem 
does not disappear in the low-lying bound states as the coupling decreases. It all returns 
to the issue of how rotational symmetry is recovered in this approach - an issue that 
extends beyond our simple leading order calculation. 

We will conclude this section with a list of approximations we make in the second step 
of the calculation: 

1. We do not include emission and absorption in Hq. 

2. We replace Xa, xj, = 1 — Xa hy Xa = ^ — r?, 1 — Xa = jj^ + m and Taylor expand 
i, in powers of r/. 

In energy denominators, we neglect terms that are higher than second order in 
momenta. The same approximation is made in the arguments of step functions. 
This leads to a "new " cutoff £, which carries only one power of transverse mass. 
It should be emphasized that the new cutoff arises only in the second step, that is, 
in the nonrelativistic approximation to the light-front effective Hamiltonian. Only 
in this context does it replace the original cutoff Elsewhere we have to work 
with ^ (e.g. the couphng runs with ^). It is which has to be at hadronic 
mass scales. In the self-energies we keep only the leading constant shift, which is 
independent of momenta. This is because the self-energies are already 0{g^). 

3. We introduce the third component of the "equal-time" momentum: kz = ^MabV^ 

23 



and extend the integration limits on Qz from {—Mat, +Mab) to (—00, +00). It is 
easy to show that this third component of the equal-time momenta coincides with 
the z-component of equal-time momentum in the lowest order in a. As mentioned 
earlier, it is possible to introduce a third component without nonrelativistic reduction 
[p!o[ , but it does not yield any simplifications if the masses are not equal. In the 
nonrelativistic limit they both agree to the leading order in powers of momenta. The 
main point we want to make here is that in the nonrelativistic limit the dynamics 
naturally lead to equal-time dynamics. 

4. For the leading order bound state calculation, we include only the rotationally in- 
variant moment of the confining potential (i.e. Vq given in (25) ) in Hq. 



This gives Hq which we choose for the purposes of the bound state calculation: 



Ho = 2M, 



ab 



2m 



Cpot CpaC 



+ 



vr 



VoiCv) 



(28) 



where Vo(^^) is given in (25) and E contains the finite shift produced by self-energies 
after subtracting terms needed to make the confining potential vanish at the origin: 



27r 



l + ^)log 



1 

+- 



AC 

ma 



mn 



C + ma 
rrib 



.,l.|^).o. 



nib 



C + mb 



1 mb 5 
AC + ma^ AC + mt ^ 2 



(29) 
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3.1 Schrodinger equation. 



We now want to find the mass of a bound state and its wave function tp{K^,x): 



J 2(27r)'^-\/a;(l — x) 



We use Lorentz invariant normalization for the states: 



{P'\P) = 2{2TTfP+6^{P - P' 



and the wave function is normalized to one: 



/ 



(Pk-^ dx 
2(27r)3 ' 



{K^,x)\' = l. 



The bound state satisfies: 



Ho\P)=M^\P) , 



(30) 



(31) 



where is the invariant mass of the bound state. Let the mass of the bound state be 



= (m„ + rribf + 2{ma + mb)E , 



(32) 



which defines E. 

Substituting for Hq and A^^ in equation (31), after some straightforward algebra one 
obtains a bound state equation for the wave function ip: 



Mab 



\ ^ X] -|- 1 

\ 2m dr^ J 



V'(r) = Mab 



1 Konf('Cr) 



^{f) , (33) 



where m is the reduced mass and we choose Konf (-Cr) = Vo{jCr) given in (25). We choose 
the confining potential for the leading order bound state calculation to be rotationally 
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invariant, but other choices are possible. To remind the reader of this possibihty and to 
keep the discussion as general as possible, we use TZ as the argument of the confining 
potential. 

It is convenient to use a dimensionless separation TZ = Cr that naturally arises in 
the confining piece of the potential, and to absorb — S into a definition of the eigenvalue 
E of the Schrodinger equation. When extracting the bound state mass, — S has to be 
subtracted. The bound state equation in the dimensionless separation is: 



^ i- + jCaCF f-V^conf (7?) + l^coul(^) 



V'(7^) = Eip{n) . (34) 



Multiplying both sides of the equation by Imj and introducing a dimensionless coupling 
and eigenvalue: 

c = , (35) 

e = (36) 

we obtain a Schrodinger equation, which depends only on dimensionless variables: 
(P /I. 



-ir + C f-Konf (^) + Koul(7^) 



V'(7e) = e^(7^) . (37) 



This form is advantageous for numerical study, but moreover, it is quite general - one 
obtains an equation of this form for any quark-antiquark systems and any choice of the 
confining potential in the nonrelativistic limit, regardless of the masses, providing they 
are nonzero. For different systems C, c, e would differ, but the resulting dimensionless 
Schrodinger equation will be the same. Thus in the leading order, qualitative characteris- 
tics of spectra depend only on one particular combination of the masses and the coupling. 
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3.2 Bohr analysis. 

The purpose of this analysis is to gain a qualitative understanding of the physics described 
by the Schrodinger eqn. (37). 

Since this is only a qualitative analysis, we will neglect the finite terms in the confining 
potential, but keep in mind that at small distances the confining potential vanishes. Since 
the Schrodinger equation is dimensionless, all quantities in this analysis are dimensionless 
also. The eigenvalue (energy e) is given by a sum of the kinetic energy, which in our case 
is simply p^, p being dimensionless; and the potential energy. We use the uncertainty 
principle to replace the momentum by ^. 

Let us consider the ground state: 



eo ~ / + F(7^) ~ + c(21og7^ - h 



(38) 



Now we find TZ which minimizes the energy: 




The solution is 



7^o = 



Vc2 + 8c - c 
2^ 



Similarly, we find the Z = 1 excited state, for which 



7^l = 



+ 3 • 8c - c 
2^ 



We now consider two limiting cases: when c is small and when c is large. 
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c small: In this case, in the ground state TZq and the energy eo ~ c + c log 

In the lowest I = 1 state TZi/TZi k, ^3, and the splitting in the energy between these two 
states is ei — eo ~ c log -\/3. 

Of course, we are interested in "real world" energies and not the dimensionless results. 
If we "unwrap" the dimensionless results, 

ei-eo ~ -^{Ei-Eq) 

c —a (39) 

which implies that 

Ei-Eo^ aJC. (40) 

The splitting between the ground state and the lowest lying P-state is independent of 
masses. 

Similarly one can show that what sets the scale for individual energies is JC and that 
the size of the system, TZ, depends both on £. and the reduced mass m. 

c large: In this case, TZq scales like |, the ground state energy eo ~ — ^, while 7?.i « | 
and ei ^ 



Unwrapping shows that the scale jC drops out: 



m „ 



^ qaf, (41) 
where e stands for any of the dimensionless energies under consideration, and E for any 
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of the "real" bindings: 

E aV (42) 

which depend only on the reduced mass but not on the confining scale. 
Similarly, jC drops out also of the expression for the size. 

In summary, when c is small, the spectra are determined by the logarithmic part of 
the potential; when c is large, the spectra are Coulombic. Recall that c is proportional 
to the reduced mass. Let us set aside questions on how C, or the original light-front 
A^, depends on the masses of constituents. There is a natural distinction between heavy 
quarkonia and other mesonic systems. In the case of heavy quarkonia the reduced mass 
is always proportional to the heavy mass, while in all other systems it is related to the 
light mass. This means that c is larger for the heavy quarkonia, and thus these states 
are substantially more Coulombic; while the other systems, including heavy mesons, are 
more influenced by the confining potential. 

This simple technique can be used to estimate expectation values of various operators 
in the ground state. 

3.3 Numerical results 

The numerical results are in agreement with the simple Bohr analysis. Figure 3 presents 
the numerical results for the dimensionless eigenvalue. It confirms that when c is small, 
the spectra are dominated by the log potential. As c increases, the Coulomb potential 
becomes more important, especially for the ground state. 
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0.2 0.4 0.6 0.8 1 



c 

Figure 3: Dimensionless eigenvalue e for the ground state, the lowest lying P state (1=1) and 
the lowest lying excited S state (1=0). 

Note that in the logarithmic regime (i.e., when c is smah) the ground state dimension- 
less energy is always larger or at least comparable to the splitting between the ground 
state and the / = 1 state. For heavy mesons the splitting is a few hundred MeV, which 
would imply that the binding in the ground state is also large. However, recall that we 
had to absorb a constant shift from the sch'-cnergies into the definition of the eigenvalue 
and that the self-energy may be fine tuned at low energies. 
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Figure 4: Ground state wavefunction compared to the Coulomb ground state wavefunction at 
the same couphng. (a) For c=0.1, the ground state wavefunction differs significantly from the 
Coulombic ground state wavefunction, while (b) for c=1.0 they are similar. 



Figure 4 shows a few typical wavefunctions for the ground state and lowest lying excited 
states. In figure 4(a), (b) we compare the exact solution of the Schrodinger equation (37) 
to the ground state wavefunction of the Coulomb problem at the same value of the couphng 
c. When c is small, the wavefunction is quite different from the Coulombic wavefunction. 
As c increases, it becomes closer and closer to Coulombic, in agreement with the Bohr 
analysis presented above. Wavefunctions of the lowest two excited states for c = 0.1 
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Figure 4: (c) Wavefunctions of the ground state, the lowest lying P state (1=1) and the lowest 
lying excited S state (1=0). 



which we later use to fit data, are shown in figure 4(c) together with the ground state 
wavefunction. 

Knowing the wavefunctions, one can calculate expectation values of various operators. 
For example, figure 5 shows the expectation value of TZ^. Even though 7?.^ >> 1 in the 
entire range of c, it does not mean that the nonrelativistic approximation is valid. Recall 
that r is related to the dimensionless TZ and is measured in units of JC; so the expectation 
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Figure 5: Expectation value of TZ"^ in the ground state. Dotted lines are connecting the data 
points. 



value of mr depends on the relation of C to the mass. 

We now attempt to unwrap the dimensionless results for B mesons. There is not 
enough experimental data, but moreover, the leading order results are too crude to justify 
a precise fit. Instead we use qualitative arguments and try to find a set of parameters 
that is reasonable. 

We know from the heavy quark effective theory that in heavy mesons the splittings 
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should be independent of the heavy masses |14|. We also know that the spectra are not 
Coulombic. Both requirements are satisfied if we choose c small. For example, with 
c = 0.1, which is on the border of the logarithmic regime, we find that we can fit 
the splitting between the two lowest lying doublets with m = 0.32 GeV, a = 0.35 and 
A = 0.98 GeV. At this value of parameters S ~ -0.37GeV, leading to EQ/m ~ 0.8. This 
shows (together with the expectation value of mr ~ 1.1) that the ground state is not 
nonrelativistic, which is not surprising. 

For systems consisting of lighter quarks, c would be smaller because the reduced mass 
decreases. As c decreases, the state becomes more and more sensitive to the confining 
potential at larger distances IZ. At some distances we expect the potential to become 
stronger than logarithmic, but without a calculation to a higher order in g we cannot 
decide whether those distances will manifest themselves in any spectra. Also, as one deals 
with lighter systems, the question of restoration of rotational symmetry becomes crucial. 



4 Summary and conclusion. 

Starting with the canonical light-front Hamiltonian with no zero modes, we use a similarity 
transformation to find an effective Hamiltonian which is band-diagonal with respect to 
a hadronic energy scale. We calculate the effective Hamiltonian to order for qq color 
singlet states with massive quarks of arbitrary masses. Then we split the Hamiltonian to 
Hq, treated to all orders, and included perturbatively, choosing the spin-independent 
and rotationally invariant part of its nonrelativistic reduction for Hq. In the nonrelativistic 
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limit the bound state equation leads to a dimensionless Schrodinger equation. Its scaling 
provides a powerful tool to classify the operators and estimate their expectation values. 

We solve the leading order problem, and find that our calculation is acceptable for 
B mesons, which we can fit with a set of reasonable, self-consistent parameters. We 
show that heavy mesons are qualitatively different from heavy quarkonia, but there are 
similarities with other mesonic systems: strange mesons and isospin 1 light mesons. Our 
approach enables us to relate different mesonic systems and use qualitative features of 
spectra (e.g., almost constant mass squared splitting in the lowest lying pseudoscalar and 
vector states, ordering of the lowest lying 0+ and 2+ states) as a check of our effective 
operators. In this paper we present the leading order problem, study of spin and angular 
momentum-dependent operators generated to this order of the coupling by the similarity 
transformation will follow. We hope that results for heavy mesons can be generalized to 
lighter mesons, at least qualitatively; although it is not clear due to violation of rotational 
symmetry. 

We manage to postpone the problem of lack of manifest rotational invariance in the 
confining potential to an order of bound state perturbation theory where there are addi- 
tional corrections beyond the current calculation. Nevertheless, the corrections are not 
suppressed by powers of the coupling and this raises a serious warning that we do not 
yet know what is needed to recover rotational invariance. There are several options - 
one might be that in our second order light-front calculation the qq approximation is 
insufficient. Since these corrections enter at the same order of bound state perturbation 
theory as the emission and absorption processes, it is possible that rotational invariance 
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requires that a qqg component be included to compensate for a rotationally noninvariant 
qq component. Another option is that this is not reahy an issue in the "real world" - the 
states which would be mixed are well separated in energy. In this case the presence of 
a rotationally noninvariant long range interactions would be merely a nuisance because 
they make it difficult to construct a simple scaling analysis that reveals the qualitative 
features of meson spectra. The resolution of these and other issues must await further 
calculations. 
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